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Abstract. The point X-ray source IE 161348-5055 is observed to display pulsations with the period 
6.67 hr and \P\ < 1.6 x 10~ 9 ss~'. It is associated with the supernova remnant RCW 103 and is 
widely believed to be a ~ 2000 yr old neutron star. Observations give no evidence for the star to 
be a member of a binary system. Nevertheless, it resembles an accretion-powered pulsar with the 
magnetospheric radius ~ 3000km and the mass-accretion rate ~ 10 14 gs _1 . This situation could be 
described in terms of accretion from a (residual) fossil disk established from the material falling 
back towards the star after its birth. However, current fall-back accretion scenarios encounter major 
difficulties explaining an extremely long spin period of the young neutron star. We show that the 
problems can be avoided if the accreting material is magnetized. The star in this case is surrounded 
by a fossil magnetic slab in which the material is confined by the magnetic field of the accretion flow 
itself. We find that the surface magnetic field of the neutron star within this scenario is ~ 10 12 Gand 
that a presence of ^ 10~ 7 Mq magnetic slab would be sufficient to explain the origin and current 
state of the pulsar. 
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1. INTRODUCTION 

The point X-ray source IE 161348-5055 (hereafter 1E1613) is observed to display 



pulsations with the period P* ~ 6.67 hr fl 1 3L 11511 . It is located near the center of the 



supernova remnant RCW 103 [|39il of the age t* ~ 2000 yr [|3jl [it] situated at the distance 
of 3.3 kpc [@, [HI]. The X-ray luminosity of the pulsar varies in the interval Lx — 
10 33 — 10 35 erg s _1 on a timescale of a few years. The X-ray emission comes from a local 
(a v ~ 600 m) region heated up to a temperature kT ~ 0.6 — 0.8 keV [|T3, UsL [l3L [l5b . The 
upper limit on derivative of the period of pulsations |P| < 1.6 x 10 _9 ss _1 has recently 
been reported by Esposito et al. [fl5ll . 

It is widely adopted that 1E1613 is a young (~ 2000 yr) neutron star. In X-rays it 
resembles an accretion-powered pulsar [11811 which accretes material at the rate [see e.g. 



25] 



071* ~ 5 x 10 13 m- 1 R 6 L 34 gs" 1 (1) 
and is surrounded by the magnetosphere of the radius 

r "^ 3xl ° 8 ^(6^)" Cm ' (2) 



Here R(, = R ns / 10 6 cm and m = M ns / 1 AMq are the radius and mass of the neutron star, 
and L34 is the average X-ray luminosity of the pulsar in units of 10 34 erg s _1 



Observations give no evidence for 1E1613 to be a close binary system Il40ll . Nev 



ertheless, Li [26] have argued that the source of the accreting material can be a fossil 
(residual) disk formed by the supernova ejecta fall-back. The lifetime of a fall-back disk 
is about 10 4 — 10 5 yr and during this time it can supply enough material to power the 
observed luminosity ifioll . 

While the conventional accretion scenario provides us with the simplest interpretation 
of the emission spectrum, it encounters major difficulties explaining the pulsar spin 
evolution. The observed spin-down behavior indicates that the spin period of the pulsar 
is smaller than (or comparable to) a so called equilibrium period, P eq , which is defined 
by equating the spin-up and spin-down torques applied to the star from the accreting 
material. If 1E1613 were a regular neutron star accreting from a Keplerian disk its 
equilibrium period would be as short as [QI |20|] 

< d) ~ 30 ^Wl-^m-^s. (3) 

Here /I30 is the dipole magnetic moment of the neutron star in units of 10 30 Gcm 3 and 
9)ti4 = Qjt*/10 14 gs _1 . Thus, to describe the pulsar spin evolution one has to assume 
that either its surface field in the current epoch is in excess of 5 x 10 15 G or the accretion 
picture differs from the Keplerian disk. The first possibility has been already discussed 
by de Luca et al. [13], Li [26] and Pizzolato et al. [1330. Here we address analysis of the 
second possibility. We show that the star under certain conditions can be surrounded by 
a fossil magnetic slab (Sect.©. Using the spin-down torque applied to the star from the 
slab (Sect. [3]) we find that the current state of the pulsar (Sect. HI) as well as its previous 
spin evolution (Sect. [5]) can be explained provided the surface magnetic field of the star 
is ~ 10 12 G. Our results are briefly discussed in Sect. [6] 



2. GEOMETRY OF THE FALL-BACK ACCRETION FLOW 



We consider a fall-back accretion [28, 42] onto a magnetized neutron star. According to 
this scenario the star of the mass M ns and the dipole magnetic moment /! after its birth 
is embedded into a dense gaseous medium of the density p^. As the star moves though 
the gas with a relative velocity v re i it captures material at a rate 971 = Ttr^p^v^i, where 
ro = 2GM ns /v 2 el is the Bondi radius. 

The material inside Bondi radius initially follows ballistic trajectories moving towards 
the star with the free-fall velocity, Vff(r) = (2GM ns /r) 1 / 2 . If it is non-magnetized and 
does not possess an angular momentum the accretion occurs in a spherically symmet- 
rical fashion on a dynamical (free-fall) timescale, % = r/vff. The flow in this case is 
decelerated by the stellar magnetic field at the Alfven radius, r& = (/l 2 /Tl\/2GM ns ) 2 ^ , 
penetrates into the magnetosphere and reaches the stellar surface flowing along the field 



lines 1240. 

A different geometry of the accretion flow can be expected if the captured material is 
frozen into a large-scale magnetic field, Bf. In this case the magnetic pressure, £ m (r) = 



B 2 (r)/Sn «r 4 , in the spherical flow increases rapidly and reaches the ram pressure, 



£ram(^) = P( r ) v ff( r ) 06 r 5 ^ 2 > at a distance I136h 

7? sh = /3 ( 7 2/3 r G (c /v rel ) 4 / 3 , (4) 

which is referred to as the Shvartsman radius [23J. Here p is the density and cq = c s (/g) 
is the sound speed in the accreting material. /3o = j8(ro), where /3 = £ t h/£m and £ t h = 
pel is the thermal pressure in the accretion flow. Studies [lH |6[ [l9|] have shown that the 
magnetized spherical flow is decelerated at the distance R s h by its own magnetic field. 
The material inside Shvartsman radius tends to accumulate in a non-Keplerian slab in 
which it is confined by the magnetic field of the flow itself. The accretion in the slab 
proceeds on the magnetic reconnection timescale, t KC = r/ r/ m VA, which under conditions 
of interest (i.e. va < Vff and r\ m <C 1) significantly exceeds the dynamical time. Here 
va = Bf/ \[A%p is the Alfven velocity in the accreting material and < r\ m < 0. 1 is the 



magnetic reconnection efficiency [37, and references therein]. The outer radius of the 
slab is limited to r ou t < ^sh and its half-thickness can be approximated by the height of 
homogeneous atmosphere, h s (r) = k# T(r) r 2 /m p GM*, where m p is the proton mass, &b 
is the Boltzmann constant and T is the gas temperature in the slab. The angular velocity 
of material in the slab in general case is limited as < Q. s \ < £1^ — (GM*//- 3 ) 1 / 2 and the 
magnetic field scales with radius as Bf (r) oc r~ 5 / 4 [f or discussion see e.g. p]. 

A formation of the magnetic slab around a neutron star accreting material from a 
magnetized wind can be expected if R s ^ > r^, which implies v re i < v ma , where 

v ma ~ 380 P l/5 m^^- Q 6/35 m% 35 cf kms- 1 (5) 



and C6 = co/ 10 cms - [|23Ll22|l. Otherwise, the magnetic field of the accretion flow can 
be neglected. 

Finally, for a Keplerian disk to form the circularization radius, r c j rc = £ 2 Q.q TqJ GM ns , 
should satisfy the condition r circ > max{r A ,i? s h}- Here £ is a dimensionless parameter 



accounting for dissipation of angular momentum in the spherical accretion flow H35H 
and = Cljrp) is the angular velocity of the material at the Bondi radius. As shown by 
Chevalier [11], the angular momentum of the material which is initially located close 
to the neutron star is insignificant. However, the outer mantel material may have a 
significant angular momentum due to its turbulization by the reverse shock wave. As 
the turbulence is subsonic the velocity of turbulent motions is v t = £cq, where < £ < 1 
and hence, = £cq/tq. The condition r c i rc > in this case can be expressed as 
v re i < v cr , where 

v„ = 4 x 10* (_|_) »/' 5 cms -. (6) 



and (^0.2 = <d/0.2 is normalized according to Ruffert [35]. 

Thus, the geometry of the fall-back accretion process onto a newly formed neutron 
star strongly depends on the physical conditions in its environment. It can be approxi- 
mated by a spherically symmetrical flow if v re i > v ma , by a Keplerian disk if v re i < v cr 
and by a magnetic slab if v cr < v re i < v ma . The latter case is discussed in this paper. 



3. SPIN-DOWN TORQUE 



The spin-down torque applied to a neutron star from the magnetic slab can be evaluated 
as [Eh 

.(si), 



\Kd I = 4nr m h s (r m ) v m p (r m ) . (7) 

Here v m = k m r m v\{r m ) is the magnetic viscosity, r m is the radius of the magne- 
tosphere of the neutron star and < k m < 1 is the efficiency parameter. v^(r m ) = 
r m [Oh — ^sl( r m)] 5 where (0 S = 2n/P s is the angular velocity of the neutron star. Finally, 
po = /i 2 nip/Ink^ T(r m ) is the gas density at the magnetospheric boundary, which is 
defined by equating the gas pressure in the slab with the magnetic pressure due to dipole 
magnetic field of the neutron star. 

Combining these parameters and taking into account that the Alfven velocity at the 
inner radius of the slab is equal to the free-fall velocity Q6|] one finds 

I (si), k m (i 2 f Q s i(r m ) ^ 

where r cor = (GMns/ft) 2 ) 1 / 3 is the corotation radius of the star. Eq. (0 represents a 
generalized form of the spin-down torque applied to a neutron star from the accreting 
material. It shows that the absolute value of the torque is limited to its conventional 
value, < \Wlco s rj^\, for r m > but can be significantly larger if the accreting material 
approaches the neutron star to a smaller distance than r^- 

The magnetospheric radius in general case can be evaluated by solving the system 

~ P{r m )c 2 s {r m ), 



LxRns 



(9) 



9ttin(r m ) 



GM n 



where 

OTin(r m ) = 4nr 5 J 4 Po (2GM ns ) l / 4 D^(r m ) (10) 

is the rate of plasma diffusion from the slab into the stellar field at the magnetospheric 
boundary and D e ff is the effective diffusion coefficient 0, The first equation in 
system © shows that the pressure of the accreting material at the magnetospheric 
boundary is equal to the magnetic pressure due to dipole field of the neutron star. The 
second equation is the continuity equation. It shows that the rate of plasma diffusion into 
the stellar field at the magnetospheric boundary is equal to the mass accretion rate onto 
the stellar surface. Solving system © and taking into account that the effective diffusion 
coefficient is limited to D e ff > DB(r m ) | fl4|l . where 

, s ck B T(r m ) 
L>B{r m ) = , — r (11) 
\6eB(r m ) 



is the Bohm diffusion coefficient, one finds that the magneto spheric radius of the neutron 
star satisfies the condition r m > r ma , where 

_ ( cml \ 2/13 ^6/13(^)5/13 

Here B(r m ) is the magnetic field strength at the magneto spheric boundary and Tq is the 
gas temperature at the inner radius of the slab. The situation r m = r ma is realized if in- 
terchange instabilities of the magneto spheric boundary are suppressed and the accreting 
material enters the stellar magnetic field due to magnetic reconnection, i.e. in a similar 
way to the solar wind penetrating into the Earth's magnetosphere [|32j]. The spin-down 
torque applied to the neutron star in this case reaches its maximum possible value (see 
EqM- 



4. MAGNETIC ACCRETION IN 1E1613 

Let us consider a situation in which 1E1613 is an isolated neutron star accreting material 
at the rate from a fossil magnetic slab. The surface magnetic field of the star in this 
case can be evaluated as 5* > 5 mm , where 

/ \ 13/6 

B mm - 7 x 10 11 T^m- 5 /% 7 / 3 4 3 ( *L ) G (13) 

V r mag / 

is the solution of Eq. (1121) for r ma = r maa andr 6 = r /10 6 K. 
The spin derivative of the pulsar, 

p2|^(sl)i 

|D|_" l^sd I 



11 2%l ' 
in this situation (i.e. r m = r mag and P s = P*) is 



3.3 x 1(T 7 s s" 1 x k m ju| IZ 1 m- 1 ' 2 (14) 



x 



45 

-3/2 



Ps\( r mag \ / Q sl (r n 
Pj \ 3x 10 8 cm/ I c& 



where Z45 = 7/10 45 gcm 2 is the moment of inertia of the neutron star. The condition 
\P*\ < 1.6 x 10 _9 ss _1 implies 



^ f2 s i(r ma ) 



Ok 



<0.005£ m \ (15) 



which may indicate that the neutron star in the current epoch rotates at almost its 
maximum possible period P mSLX = 2n / £l s \(r ma ) 0]. 



5. POSSIBLE HISTORY 



The spin-down time, T ~ P S /2P, of a neutron star accreting material from the magnetic 
slab in the general case can be expressed as 

' ~ — , (16) 



2P A 2\x 

where 

|^( sl )|p2 p 2 k „i 

I sd rs _ . ' s K mH 



p , = 1 sq 1 ° ~ ° " ,r fn^ 



is the spin-down rate evaluated under the assumption £l s \ <C co s . Putting parameters of 
1E1613 derived in the previous section to Eq. ([TBI one finds 

r a * 1880 M3o 17/1 ^ 8/13 ^r 6 3/13 3jl- 6/13 yr. (18) 

This shows that a neutron star with the surface magnetic field of ~ 10 12 G would have 
enough time to slow down to a long period if it accretes material from the magnetic slab 
at an average rate ~ 10 14 gs _1 . The spin period which the neutron star is able to reach 
during this time is limited to P max < 2n/Q. s \(r m ) and, therefore, is determined by the 
angular velocity of the material at the inner radius of the slab. 

Studies of pulsar population [see e.g., and references therein] suggest that 

the average birth spin period of neutron stars is Pq ~ 0.5 s. The ejector (spin-powered 
pulsar) state under these conditions can be avoided if the initial mass-transfer rate in the 
fall-back accretion process is SOTo > 9jt e j, where 

9)t ej ^8xl0 13 / m M3W (J^J gs- 1 (19) 
is the solution of the equation 7^ = P e j, and 

7> ej ~ 0.26 /i /4 /i^ 2 9ft7 5 1/4 Vg 1/4 s (20) 

is the spin period at which the pressure of relativistic wind ejected from the magneto- 
sphere of a newly formed neutron star is equal to the ram pressure of the surrounding 



material at the Bondi radius [for discussion see e.g., |21j]. Here vs = v re i/10 8 cms 1 



DJl 15 = QHq / 10 15 g s 1 and f m = 1 + sin 2 y, where % I s me angle between the magnetic 



and rotational axes of the neutron star 038]]. 

A higher value of WIq would be required for a newly formed neutron star to avoid the 
propeller state and to start its spin evolution as an accretor. This situation can be realized 
if Otto > SEKpr, where 

/ d \ -13/6 

M^sxm'Vi'V/^A;) g s-' (2i) 



is the solution of the equation Pq = P pr , and 



P pr ~ 3.5 >n- 5/l3 T 6 3/U m H 6/n s (22) 



is the spin period defined by equating r ma — r cor [see Eq. 19in|2lH. 

Thus, one can envisage the following evolution scenario of 1E1613. The neutron star 
was formed in the core-collapsed supernova explosion with a surface field Bq ~ 10 12 G 
and a spin period Pq ~ 0.5 s. It started its spin evolution in the accretor state being 
surrounded by a magnetic slab formed by the supernova ejecta fall-back. The mass- 
transfer rate in the slab initially was in excess of 3 x 10 15 gs _1 and had been decreased 
over 2000 yr by almost two orders of magnitude. The spin period of the star during 
this time has been increased to its presently observed value defined by the angular 
velocity of the material in magnetic slab at the magneto spheric radius. The star in current 
epoch undergoes accretion from the magnetic slab at the rate and the radius of its 
magnetosphere is r ma ~ r#. The surface field of the star in the current epoch is about 
its initial value. The spin-down torque applied to the star from the accreting material is 
given by Eq. (O in which ^ s i(r ma ) ~ 0.995 k~ l co s . 



6. DISCUSSION 

In the frame of the magnetic accretion scenario 1E1613 appears to be a regular isolated 
neutron star with the conventional values of basic parameters. The extremely long period 
of the pulsar in this approach is associated with peculiar physical conditions in the stellar 
environment. Namely, it is assumed that the neutron star is surrounded by a (residual) 
fossil magnetic slab. This implies that the supernova ejecta was magnetized and the 
relative velocity of the material moving towards the neutron star after its birth met the 
condition v cr < v re i < v ma (see Eqs.[5]and[6]). 

Our basic assumption seems rather plausible in the light of current views on the fall- 
back accretion process. Numerical simulations of supernova explosion [see e.g.. l29l[l6|l 
favor a strong magnetization of the ejecta in the vicinity of a newly formed neutron 
star. Our assumption is also consistent with the results of studies of supernova remnants 
[see e.g. 0, and references therein] which suggest that the magnetic pressure in the 
environment of the neutron stars is comparable to the thermal gas pressure, i.e. /3 ~ 1. 
The velocity of the material falling-back towards the neutron star after its birth lies in 
the range 200 — 2000km s _1 lf4lL fllTl . This is substantially larger than the typical value 
of v cr (see Eq. [6]), but is comparable with v ma provided the temperature in the material 
in the vicinity of the star after its birth is ~ 10 6 — 10 7 K. 

Finally, the initial mass-transfer rate in the slab, 9JT ~ 10 15 — 10 16 gs _1 , evaluated 
in our scenario meets the conditions previously reported by Chatterjee, Hernquist & 
Narayan [10]. Following their results one can expect that the mass-transfer rate in the 
slab has been decreased over a time span of ~ 2000 yr by almost 2 orders of magnitude 
[see also 0] to its current value 10 13 — 10 14 gs _1 . This indicates that a presence of a 
magnetic slab of the mass M s \ Wl pr T* ~ 10~ 7 Mq would be sufficient to supply enough 
material for powering the X-ray luminosity of the pulsar over 2000 yr by the accretion 
process. 
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